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BACKGROUND 

For banks and other financial institutions, risk measurement plays a central role. Risk 

levels must conform to the capital adequacy rule. An error in the computed risk level may 

thus affect a bank's investment strategy. The state of the art is measuring risk by analyzing 
15 daily data: using one market price per working day and per financial instrument. In this 

description, the stochastic error of such a risk measure is demonstrated in a new way, 

concluding that using only daily data is insufficient. 

The challenge for statisticians is to analyze the limitations of risk measures based on 

daily data and to develop better methods based on high-frequency data. This description 
20 meets this challenge by introducing the time series operator method, applying it to risk 

measurement and showing its superiority when compared to a traditional method based on 

daily data. 

Intra-day, high frequency data is available from many financial markets nowadays. 
Many time series can be obtained at tick-by-tick frequency, including every quote or 

25 transaction price of the market. These time series are inhomogeneous because market ticks 
arrive at random times. Irregularly spaced series are called inhomogeneous, regularly spaced 
series are homogeneous. An example of a homogeneous time series is a series of daily data, 
where the data points are separated by one day (on a business time scale which omits the 
weekends and holidays). 

30 Inhomogeneous time series by themselves are conceptually simple; the difficulty lies 

in efficiently extracting and computing information from them. In most standard books on 
time series analysis, the field of time series is restricted to homogeneous time series already 
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in the introduction (see, e.g., Granger C.W.J, and Newbold P., 1977, Forecasting economic 
time series. Academic Press, London; Priestley M.B., 1989, Non-linear and non-stationary 
time series analysis. Academic Press, London; Hamilton J.D., 1994, Time Series Analysis, 
Princeton University Press, Princeton, New Jersey) (hereinafter, respectively. Granger and 
5 Newbold, 1977; Priestley, 1989; Hamihon, 1994). This restriction induces numerous 
simplifications, both conceptually and computationally, and was almost inevitable before 
cheap computers and high-frequency time series were available. 

SUMMARY 

10 U.S. Provisional Application No. 60/200,743, filed May 1, 2000, discloses a new time 

series operator technique, together with a computationally efficient toolbox, to directly 
analyze and model inhomogeneous as well as homogeneous time series. This method has 
many applications, among them volatility or Value-at-Risk (VaR) computations tick by tick. 
A comparison is made herein between VaR results based on daily data, sampled at a 

15 certain daytime, and results based on tick-by-tick data and the new time series operator 

technique. If using daily data, a surprising and (for practitioners) alarming sensitivity against 
the choice of the sampling daytime is observed. The stochastic noise seems higher than 
acceptable to risk managers. An alternative VaR computation based on tick-by-tick data and 
a new time series operator technique is shown to have similar properties, except for two 

20 advantages: distinctly reduced noise and availability of up-to-date results at each tick. 

The time series operators can also be used in the formulation of old and new 
generating processes of time series. This opens new ways to develop process equations with 
new properties, also for inhomogeneous time series. 

A preferred embodiment comprises a method for determining value-at-risk based on 

25 tick-by-tick financial data. Major steps of the method comprise the following: (1) financial 
market transaction data is electronically received by a computer; (2) the received financial 
market transaction data is electronically; (3) a time series z is constructed that models the 
received financial market transaction data; (4) an exponential moving average operator is 
constructed; (5) an operator is constructed that is based on the exponential moving average 

30 operator; (6) a causal operator Q.[z\ is constructed that is based on the iterated exponential 
moving average operator; (7) values of predictive factors are calculated; (8) the values 
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calculated by the computer are stored in a computer readable medium, and (9) value-at-risk is 
calculated from the values stored in step (8). 

BRIEF DESCRIPTION OF THE DRAWINGS 
5 FIG. 1 illustrates the relationship between a time series Z and a time series operator 

Q. 

FIG. 2 depicts an example of a causal kernel co(0 of a moving average. 
FIG. 3 depicts a graph of a kernel of a simple EMA operator. 
FIG. 4 depicts graphs of selected EMA operator kernels. 
10 FIG. 5 depicts graphs of selected MA operator kernels. 

FIG. 6 depicts graphs of selected terms of a kernel of a differential operator A. 
FIG. 7 illustrates volatility of standard RiskMetrics. 
FIG. 8 illustrates operator-based tick-by-tick volatility. 
FIG. 9 compares RiskMetrics to operator-based volatility. 

15 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
1. The time series operator technique 

In this description, only a minimum of a description of time series operators is given, 
so the applications of the following sections can be understood. The theory of the time series 
, 20 operators is explained in U.S. Provisional Application No. 60/200,743 
1. 1 Inhomogeneous time series 

A time series z consists of elements or ticks z, at times t,. The sequence of these time 
points is required to be growing, t,> t,_-^. 

A general time series is inhomogeneous, meaning that the sampling times 
25 are irregular. For a homogeneous time series, the sampling times are regularly spaced, - 
= St. 

For some discussions and derivations, a continuous-time version of z has to be 
assumed: z(f). However, the operator methods that are eventually applied only need the 
discrete time series (t, , z,). 
30 The letter x is used to represent the time series of logarithmic middle prices, x = 

(In /7yd AskV^. This quantity is used in the applications. 
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1.2 Operators 

An operator Q., mapping from the space of time series into itself, is depicted in FIG. 
1. The resulting time series n[z] has a value of iQ[z](?) at time t. Important examples are 
moving average operators and more complex operators that construct a time series of 
5 voliatility from a time series of prices. 

Linear and translation-invariant operators are equivalent to a convolution with a 
kernel a> : 

Cl[z]it}=l'Giit-t')z(t)dt' (1) 

10 A causal kernel has co(0 = 0 for all / < 0. No information from the "future" is used. If 

G)(/) is non-negative, Q.[z] is a weighted moving average of z whose kernel should be 
noimalized: 

J'co(/-/Odr'=|^°°co(0dr = 1 (2) 

15 This kernel co(0 is the weighting function of the past. FIG. 2 depicts an example of a causal 
kernel &{t) of a moving average. 

The range of an operator is the first moment of its kernel: 

r = j"° co(t)tdt (3) 

20 Thiis indicates the characteristic depth of the past covered by the kernel. 

Operators are useful for several reasons, as will be shown. One important aspect is to 
replace individual ticks from the market by local short-term averages of ticks. This mirrors 
the view of traders who consider not only the most recent tick but also the prices offered by 
other market makers within a short time interval. 
25 1.3 The simple EA£4 operator 

The Exponential Moving Average (EMA) operator is a simple example of an 
opeirator. It is written EMA [r; z] and has an exponentially decaying kernel (as shown in FIG. 
3 vs^hich depicts a graph of the kernel, with -r=l): 

co(/) = ema(t) = 

30 T 

(4) 
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According to eqs. (3) and (4), the range of the operator EMA[t; z] and its kernel is 

(5) 

5 Th; variable z thus characterizes the depth of the past of the EMA. 

The values of EMA[t; z](t) can be computed by the convolution of eq. (1), if z(t) is 
known in continuous time. This implies an integration whose numerical computation for 
many time points t is costly. Fortunately, there is an iteration formula that makes this 
computation much more efficient and, at the same time, solves the problem of discrete data. 
1 0 Thi s means that we do not need to know the whole function z(t); we just need the discrete 
time series values z, = z(r,) at irregularly spaced time points /,. The EMAs are calculated by 
the following iteration formula: 

ElV[A[r; z](0 = ^ EMA[r; z] (t„_,) + (I - ju) z (t„) + (m - v)[z(t„) - z(t„_,)] (6) 
15 with 

L-fu 1 

M=e~^ a=^^-^ (7) 

T 

(8) 

a 

20 Thi s variable v is related to the problem of using discrete data in a convolution defined in 
coEitinuous time. We need an assumption on the behavior of z(t) between the discrete time 
points t,. Eq. (8) is based on the assumption of linear interpolation between points; other 
formulas for v are implied by other interpolation assumptions, as explained in U.S. 
Provisional Application No. 60/200,743. In the case of assuming the value of the old tick for 

25 the whole interval before the new tick, the correct formula is 

v=l (9) 

For a homogeneous time series, ju and v are constants. A homogeneous time series 
can alternatively be regarded as a truly discrete time series to which interpolation does not 
30 apply. This is mentioned here because it is a popular approach used by traders. For such a 
dis(jrete time series, — t„_j is defined to be 1, and the following definition is appropriate: 
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H ^ V = — — = = — ^ (10) 

1 + a T + /'„-r„_i r + 1 

The range of an operator for a genuine discrete time series has a new definition: 
r 5^7=0 ^! For EMA, this means r = ju/(l - fi) = x with OO; = (1 - The // and v 
5 values resulting from eq. (10) are very similar to those of eqs. (7) and (8) as long as a is 
small. 

The iteration equation (6) is computationally efficient, extremely so when compared 
to £L numerical convolution based on eq. (1). No other operator can be computed as 
efficiently as the simple EMA operator. However, there are means to use the iteration 
10 equation (6) as a tool to efficiently compute operators with other kernels, as shown below. 

An iteration formula is not enough. We have to initialize EMA [r; z] (?o) at the start of 
the time series. For this, we can take Zq = z (to) or another typical value of z. This choice 
introduces an initial error of EMA [t; z] (Q which decreases exponentially with time. 
Therefore, we also need a build-up period for EMA [t; z]: a time period over which the 
1 5 vahies of EMA [t; z] should not yet be applied because of their initial error. Build-up periods 
should be multiples of t, e.g., 5r. The choice of a large enough build-up period is discussed in 
U.S. Provisional Application No. 60/200,743. 

1.4 The operator EMA [t, n; z] 

Time series operators can be convoluted: a time series resulting from a an operator 
20 can be mapped by another operator. This is a powerful method to generate new operators 
witii different kernels. 

The EMA [t, n; z] operator results from the repeated application of the same simple 
EMA operator. The following recursive definition applies: 

EMA {% k; z] = EMA [r ; EMA {x, 1; z]] (11) 

25 

with EMA [r, 1; z] = EMA [r; z]. The computationally efficient iteration formula of the 
simple EMA, eq. (6), can again be used; we have to apply it recursively {n times) for each 
n^w tick (tp Zj). For v, we insert eq. (8) which is based on a linear interpolation assumption 
betAveen ticks. (This assumption is just a good approximation in some cases, as discussed in 
30 U.S. Provisional Application No. 60/200,743. 

The operator EMA [t, n; z] has the following kernel: 
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ema[r,«](0 = M- ^ (12) 

(n - 1)! t) t 

Thiis kernel is plotted in FIG. 4, which depicts graphs of selected EMA operator kernels, for 

se-veral n(n = \ (thin), 2, 3, 4, and 1 0 (bold); left graph is for t = 1 , right graph is for r = « t = 
5 1). For large n (e.g., « = 10 in FIG. 4), the mass of the kernel is concentrated in a relatively 

narrow region around a time lag of nr. The corresponding operator can thus be seen as a 

smoothed backshift operator. 

The family of functions of eq. 12 is related to Laguerre polynomials which are 

ortliogonal with respect to the measure e' (for t = 1). 
10 Operators, i.e., their kernels, can be linearly combined. This is a powerful method to 

generate more operators. Linear combinations of EMA[t, n; z] operators with different n but 

identical t values have kernels that correspond to expansions in Laguerre polynomials. This 

means that any kernel can be expressed as such a linear combination. The convergence, 

however, of the Laguerre expansion may be slow. 
15 In practice, a small set of useful operators can be prepared with all the kernels needed. 

Asiide from the discussed expansion, it is also possible to linearly combine kernels with 

different z values. Some useful types of combined operators are presented in U.S. Provisional 

AppHcation No. 60/200,743. 

1.5 The operator MA[t, n; z] 
20 The moving average (MA) operator has kernels with useful properties as shown in 

FIG. 5, which depicts graphs of selected MA operator kernels (w = 1, 2, 4, 8, and 16; t = 1). 

It is constructed as a sum of EMA [t, n; z] operators: 

UA[t,n] = - VEMA[t',A:] with r' = — (13) 

n k=i n+\ 

25 

Thi3 variable r 'is chosen such that the range of MA[t, «] is r = r, independent of n. For « = 1, 
we obtain a simple EMA operator, for « = °° the rectangularly shaped kernel of a simple 
moving average with constant weight up to a limit of It. This simple rectangular moving 
avc;rage has a serious disadvantage in its dynamic behavior: additional noise when old 
30 observations are abruptly dismissed from the rectangular kernel area. Kernels with finite n 
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are better because of their smoothness; the memory of old observations fades gradually rather 
than abruptly. 

The formula for the MA[t, n] kernel is 



«(o.iiii£-!ii:±[x]' (14) 

Many other kernel forms can be constructed through different linear combinations of 
E]V[A[t, n; z] and other operators. 

1. 6 From returns to differentials 

Most statistics in finance is based on returns: price changes rather than prices. 
Simple returns have a rather noisy behavior over time; we often want differences between 
local averages of x: smoothed returns. 

Smoothed returns are computed by differential operators. Examples: 

X - EMA[t, n; x], where the EM A replaces x (t-r). This is used by the 
application of section 3.2. 

EMA[r,, «,] - EMA[t2, WjI^ with t, < Tj or < «2 • 

A[t] = y {EMA[aT, 1 ] + EMA[aT, 2] - 2 EMA[«;eT, 4] } , with y = 1 .22208, p = 
0.65 and a' = y(8p - 3). The is normalized, so A[t, 1] = 0, A[t, t] = 1. The 
kernel of this differential operator, described in U.S. Provisional Application 
No. 60/200,743, is plotted in FIG. 6, which depicts graphs of selected terms 
of a kernel of a differential operator A (full curve is kernel of A[t] (t = 1); 
dotted curve is the first two terms of the operator; dashed curve is the last term 
2 EMA[a/3r, 4] ) . 

The expectation value of squared smoothed returns may differ from that of the 
corresponding simple returns. This has to be accounted for when comparing the two 
concepts, for example in terms of a factor c in eq. (20). 

1. 7 Volatility measured by operators 

Volatility is a central term in risk measurement and finance in general, but there is no 
unique, universally-accepted definition. There are volatilities derived from option market 



prices and volatilities computed from diverse model assumptions. In this description, the 
focus is on historical volatility: a volatility computed from recent data of the underlying 
instrument with a minimum of parameters and model assumptions. 

For computing the time series of such a volatility, a time series operator is again the 
suitable tool. We first define the nonlinear moving norm operator: 

MNorm[T,;7;z] = {MA[T;izr]}''^ (15) 
Thiis operator is based on a linear MA operator (w^here we are free to choose any positive, 
causal kernel); it is nonlinear only because a nonlinear fimction of the basic time series 
variable z is used. MNorm[T, p; z] is homogeneous of degree 1 . 

The volatility of a time series x can now be computed with the help of the moving 
noim operator: 

Volatility [t^,T2,p;x] = MNorm 



— ,/?;A[t2;x] 



(16) 
(17) 



This is the moving norm of (smoothed) returns. With = 2, it is a particular version of the 
frequently used RMS value. However, some researchers had and have good reasons to 
choose a lower value such as = 1 in their special studies. 

Eq. (17) is based on a moving average (MA) and a differential (A) operator. In 
principle, we may choose any MA and A operator according to our preference. In the 
applications of section 3, this choice is made explicit. 

The volatility definition of eq. (17), as any definition of historical volatility, 
necessarily has two timing parameters: 

1 . the size of the return measurement intervals: x^'^ 

2. the size of the total moving sample: Tj, often » Xj; defined as the double range 
of the used MA. The MA operator has a range (center of gravity of the kernel) 
ofTi/2. 



30 2. Application: volatility computation in risk management 
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Computing recent volatility is a central ingredient of risk assessment in risk 
management. Here it serves as an example to demonstrate the usefulness and superiority of 
time series operators. 

The RiskMetrics^"^ method (see J.P. Morgan, 1996, RiskMe tries - technical 
5 document. Technical report, J.P. Morgan and International marketing - Reuters Ltd.) is 
chc»sen as a well-known example. First it is shown to be a special application of the time 
seriies operator technique. Then a better volatility computation method, also based on the 
time series operator technique, is proposed. Thus two approaches are compared: 

1 . The RiskMetrics method, based on an IGARCH model with working-daily 
10 data. 

2. A tick-by-tick alternative, following RiskMetrics as closely as possible, based 
on time series operators. 

In both cases, squared volatility is defined as the expectation tr ^ of squared, working- 
daily changes of the logarithmic middle price x. This is for the sake of a meaningful 
1 5 comparison; it does not imply that using squared returns is necessarily the best choice for an 
optimal volatility definition. 

2. 1 Conventional computation based on daily data 

The RiskMetrics method is based on an IGARCH model. Its volatility formula gives 
the conditional expectation of the squared return assuming IGARCH: 
20 o^it) = o^(t - 1 wday) + ( 1 - //) [x{i) -x{t-\ wday)]^ ( 1 8) 

vnXh. \i = 0.94. This is just an EMA iteration which can also be written in our operator 
notation: 

(7^(0 = EMa[x = 15.67 wdays; [x{t) -x{t - 1 wday)]^] (19) 
ev£iluated at discrete time points separated by 1 working day (= 1 wday); with EMA range t = 
25 ///(I - fj.) in working days, following eq. (10). 

Thanks to the regularity of the underlying homogeneous time series, fi = 0.94 is a 
constant. In general, the constancy of // makes the operator technique particularly efficient 
for homogeneous time series. 

FIG. 7, which illustrates volatility of standard RiskMetrics, shows the resulting 
30 volatility as a function of time, in an empirical example from the foreign exchange (FX) 

market: USD/JPY data in January and February 1999. The volatility is computed only once 
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per working day, at a given time of day; the resulting volatility value is valid until it is 
replaced by a new one, one working day later. Circles show data sampled at 7 am GMT, and 
diamonds show data sampled at 5 pm GMT. Computations are independent. The price is 
plotted against time on the lower graph. 
5 In FIG. 7, two such volatilities are plotted. The difference between the two curves 

solely originates from the choice of time when the raw data x is sampled and the volatility is 
coraputed by eq. (1 8) or (19). One curve is sampled at 7 am GMT which is a time in the late 
afternoon of East Asian time zones - a suitable daytime for the daily risk calculations of an 
East Asian risk manager. The other curve is sampled at 5 pm GMT - a suitable daytime for a 
1 0 risl: manager in London. 

The differences between the two curves are surprisingly large: up to 25%, an alarming 
unc;ertainty for risk managers. Risk levels are linked to a bank's capital through the capital 
adequacy rule, so differences in risk measurements have a major impact on banking. In our 
case, two risk managers measiire very different volatility and thus risk levels for the same 
1 5 finimcial instrument, just because they live in different time zones. A difference can persist 
ove;r weeks, as shown in FIG. 7. This figure is just an example. The same surprisingly 
strong effect can be foimd also for other financial instruments, sampling periods, and choices 
of time of day for sampling. 

Both deviating volatility values cannot be right at the same time; there must be an 
20 error in these values. This error is of a stochastic nature; there is no systematic bias 

dependent on the daytime. In FIG. 7, the difference between the two curves is neither always 
positive nor negative; it changes its sign. 

FIG. 7 demonstrates the large stochastic error of the RiskMetrics method. The large 
size of this error has two main reasons: 
25 1 . The rather small range of the kernel of some 16 working days. The number of 

independent observations is limited. We caimot essentially change this fact, 
because the choice of a short range is also motivated by the goal of fast 
adaptivity to new market events. 
2. The results depend on only one observation per day, taken at a certain time. 
30 All the other information on prices of the day is thrown away. The value at 

that time may be little representative for the full day: it may be located on top 
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of a short-lived local peak of the price cxirve. this is indeed the reason for the 
large deviations of the two curves in FIG. 7. The effect is exacerbated by the 
known fact that returns have a heavy-tailed distribution function: extreme 
(intra-day) events dominate the statistics. 
5 The focus here is not so much the behavior of RiskMetrics (IGARCH), but the 

problems of using homogeneous, daily data in general, no matter which GARCH-type or 
other model is investigated. The significance of most results can be improved by using all the 
available information, tick by tick, as shown in the next section. 
2. 2 Tick-by-tick volatility computation 
1 0 For the sake of a fair comparison, a tick-by-tick volatility computation is introduced 

that follows RiskMetrics as closely as possible. There are two innovative modifications: 

The squared volatility a\t) is computed at every available tick, not just once 
per working day. 

Simple returns are replaced by operator-based, smoothed returns. 
1 5 Nothing is changed otherwise; the sampling range of 1 5.67 working days and the working- 
daily nature of (smoothed) returns are preserved. 

The new volatility measure is again defined in operator notation (where "wdays" 
stands for working days): 

= cEMA[T=15.67wdays;(x-EMA[T = lwday,4;x])2] ^20) 
20 Th is is just a special case of eq. (16). The computation is efficiently done at every new tick, 
repeatedly using the iteration formula (6). This works not only for the simple EMA but also 
for EMA[x, 4; x] as explained in section 2.4. 

The constant c compensates for the fact that we use smoothed returns x - EMA[t, 4; x\ 
as introduced in section 2.6 instead of the simple returns of section 3.1. In the case of x 
25 following a Gaussian random walk, the theoretically correct value is c = 128/93 . Using this 
factor eliminates a systematic bias of the tick-by-tick volatility as compared to the 
RiskMetrics volatility. 

Eq. (20) is computed on a special business time scale defined as follows. The 49 
weekend hours from Friday 8 pm GMT to Sunday 9 pm GMT are compressed to the 
30 equivalent of only 1 hour outside the weekend. This fully corresponds to the time scale of 
RiskMetrics which omits the weekend days. A more sophisticated and appropriate choice of 
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the business time scale would be the ^-time of Dacorogna et al. (1992) <Dacorogna M.M. 
MuUer U.A., Nagler R.J., Olsen R.B., and Pictet O.V., 1993, A geographical model for the 
daily and weekly seasonal volatility in the FX market. Journal of International Money and 
Finance, 12(4), 413-438), but this is avoided here in order to keep the approach as close to 
5 RiskMetrics as possible. 

FIG. 8, which illustrates operator-based tick-by-tick volatility, shows the resulting 
volatility as a function of time. The same financial instrument and sampling period is studied 
as in FIG. 7. High-frequency data is available here. Now, the large differences between 
values at 7 am GMT and 5 pm GMT have vanished. The observations at these times appear 

10 as points on one continuous, consistent curve. In fact, we can obtain volatility values at any 
time of day now, not just once or twice a day. A risk manager in London essentially 
measures the same risk of the instrument as a risk manager in East Asia, as should be 
expiected in normal situations. The risk levels deviate only if a dramatic event between the 
two daytimes of measurement happens. This is natural; the operator-based volatility quickly 

1 5 reacts to dramatic events, as can be seen in FIG. 8. 

The variations of the volatility level over time are moderate in FIG. 8. The extreme 
volatility minima and minima of FIG. 7 which are mostly due to stochastic noise have 
vanished. The new tick-by-tick volatility has less stochastic noise than the RiskMetrics 
volatility, although the moving sample range of 15.67 working days is the same. 

20 The curves of FIGS. 7 and 8 are combined in FIG. 9, which compares RiskMetrics to 

operator-based volatility. Here, we can see that the tick-by-tick volatility (bold curve) has a 
dynamic behavior similar to the known volatilities while avoiding extreme oscillations due to 
stochastic noise. 

The lower noise level of the tick-by-tick volatility is now plausible, but we need 
25 sciisntific evidence for this. In the general case, such evidence can be gained through Monte- 
Carlo studies based on a certain process assiimption, comparing the error variances of the 
RiskMetrics volatility and the tick-by-tick volatility. In the case of a Gaussian random walk, 
we have even stronger evidence: by using continuously-overlapping returns instead of non- 
ov(jrlapping returns, the error variance of the empirically determined cr^ is reduced to 2/3 of 
30 the original value. 
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The tick-by-tick operator is indeed using (almost) continuously overlapping returns. 
In addition to this, it is based on smoothed rather than simple returns, which also leads to a 
reduction of stochastic noise. 

Other advantages of tick-by-tick, operator-based methods are the efficient 
5 coraputation based on iterations and the updating at every new tick. Thanks to fast updates, 
the volatility measure can quickly react to new market events such as shocks, at any time of 
day. 

2. 3 VaR computation in real-time 

Conventional Value-at-Risk (VaR) computations are done once a day, usually in the 
1 0 evening. Traders and portfolio managers typically do not know their current risk; they just 
know yesterday evening's risk. What they really need is a real-time VaR computation, 
updated as quickly as they can change their positions. 

The tick-by-tick operator proposed in section 3.2 and eq. (20) is a tool to make a real- 
time VaR possible. 

15 A real-time VaR computed according to these guidelines would still be somewhat 

similar to RiskMetrics, except for the substantial benefits of lower noise and a higher 
updating frequency. There are many criticisms of RiskMetrics that would apply to it, too. 
Soime researchers, for example, replace the IGARCH model by another GARCH-type model. 
Otlier researchers focus on the behavior of extreme price changes which may follow other 

20 lav/s than average-size changes. Moreover, return observations over intervals other than daily 
(for example, hourly or weekly) returns contain valuable information that should also be used 
in a VaR computation. 

3. Process equations based on operators 

Processes of the ARM A and GARCH families can be expressed in terms of time 
25 series operators, as we have seen for IGARCH in eq. (19). The squared conditional volatility 
of a GARCH(1,1) process, for example, can be written as follows: 

cy2(0 = c + aa'it') + b[x(t)-x(t' -At)f (21) 
whiere t'=t- At. The following alternative notation is based on a simple EMA operator: 



^(t) = + EMA — ; [x(t ') - ht)f (22) 

\-a \-a \\-a I 
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This rephrased form of GARCH(1,1) is a starting point of interesting new developments. 
Initially, it applies to a discrete, homogeneous time series in the sense of eq. (10), but it 
allows for a direct and efficient computation of the GARCH(1,1) volatility from 
inhomogeneous data, since the operator technique is also suited to inhomogeneous time 
series. 

Moreover, eq. (22) can be modified to obtain other processes. The kernel of the EMA 
operator can, for example, be replaced by other kernels. The return x{t ")-x{t'- At) can be 
replaced by a smoothed return computed by a differential operator which reflects the 
perception of market participants better that the simple return. 

Dacorogna et al. (1998) (Dacorogna M.M., Miiller U.A., Olsen R.B., and Pictet O.V., 
1998, Modelling short-term volatility with GARCH and HARCH models, published in 
"Nonlinear Modelling of High Frequency Financial Time Series" edited by Christian Dunis 
and Bin Zhou, John Wiley, Chichester, 161-176) have introduced the EMA-HARCH process 
to model some empirical facts of high-frequency data in finance: the long memory of 
vol atility, the fat tails of the return distribution and the asymmetric causality between fine- 
grained (high-resolution) and coarse-grained (low-resolution) volatilities as found by Miiller 
et al. (1997) (MuUer U.A., Dacorogna M.M., Dave R.D., Olsen R.B., Pictet O.V., and von 
W(;izsacker J.E., 1997, Volatilities of different time resolutions - analyzing the dynamics of 
market components. Journal of Empirical Finance, 4(2-3), 213-239). This is one of the first 
processes whose equation is written with the help of a time series operator: 

^(0 = c.±c/jit) (23) 

The "partial volatilities" correspond to market segments and are written in terms of 
the: EMA operator: 



(5j{t) = EMA 
mihk, = 1 andA:j = 4^-^+ 1 fory> 1. 



ili— i At; [x(t ') - x(t '-kjAt)f 



(24) 



4. Conclusions 

Most financial markets produce inhomogeneous data, irregularly spaced in time. The 
time series operators described herein are able to directly use inhomogeneous data to estimate 
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statistical variables such as volatilities. This computation is made efficient by using iteration 
fonnulas. The operator technique efficiently works also for homogeneous, equally spaced 
data. 

Starting from the simple exponential moving average (EMA) operator, large families 
5 of operators with different kernels and different purposes can be constructed. A wider 
ove rview of these more complex operators, which are still computationally efficient, is 
described in U.S. Provisional Application No. 60/200,743. One example is a tick-by-tick 
Fourier analysis on a moving sample. 

Thanks to averaging, the operator technique often produces results with less noise 

1 0 (lower stochastic errors) than conventional methods based on homogeneous time series. This 
is a lso the case for the main application of the above description: volatility of daily returns as 
needed for Value-at-Risk (VaR) computations. The conventional RiskMetrics methods have 
a rather high stochastic error which is demonstrated in a new way: volatility is computed 
twice with daily data. In one case, the data is always sampled at 7 am GMT (late afternoon in 

1 5 East Asia), in the other case at 5 pm (late afternoon in London). The results of the two 

computations can differ by some 25% for many days in a row - an embarrassing fact for risk 
managers. The tick-by-tick alternative of a preferred embodiment based on time series 
operators does not have this sensitivity against the choice of the sampling time of day and has 
les s noise, while keeping the essential preferred characteristics of the known RiskMetrics 

20 method: it is still based on daily returns and still has a range (center of gravity of the kernel) 
of around 16 working days. 

The same technique is preferably used to determine a real-time VaR, updated with 
ev<jry new incoming tick from a market, with less noise than the corresponding results from 
conventional methods. Many methods of calculating VaR from volatility are known in the art 

25 (see, for example. Chapter 14: Value at Risk of Options, Futures, and Other Derivatives, by 
Jolm C. Hull (4^^ ed. 2000). 

Finally, the operator technique can be used to formulate time series generation process 
equations. This is possible for well-known processes such as GARCH and new, more 
complex processes such as EMA-HARCH. The formulation in terms of operators has many 

30 advantages: flexibility and applicability to irregularly-spaced time series. 
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Although the subject invention has been described with reference to preferred 
embodiments, numerous modifications and variations can be made that will still be v^thin the 
sco]pe of the invention. No limitation with respect to the specific embodiments disclosed 
herein is intended or should be inferred. 
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